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A  high  efficient  assembly  technique  for  large  proton  exchange  membrane  fuel  cell  (PEMFC)  stacks  is 
proposed  to  obtain  the  optimal  clamping  load.  The  stack  system  is  considered  as  a  mechanical  equivalent 
stiffness  model  consisting  of  numerous  elastic  elements  (springs)  in  either  series  or  parallel  connections. 
We  first  propose  an  equivalent  stiffness  model  for  a  single  PEM  fuel  cell,  and  then  develop  an  equivalent 
stiffness  model  for  a  large  PEMFC  stack.  Based  on  the  equivalent  stiffness  model,  we  discuss  the  effects 
of  the  structural  parameters  and  temperature  on  the  internal  stress  of  the  components  and  the  contact 
resistance  at  the  contact  interfaces,  and  show  how  to  determine  the  assembly  parameters  of  a  large  fuel 
cell  stack  using  the  equivalent  stiffness  model.  Finally,  a  three-dimensional  finite  element  analysis  (FEA) 
for  a  single  PEMFC  is  compared  with  what  the  equivalent  stiffness  model  predicts.  It  is  found  that  the 
presented  model  gives  very  good  prediction  accuracy  for  the  component  stiffness  and  the  clamping  load. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cell  technology  has  been  achieved  a  great  jumping  develop¬ 
ment  during  the  past  decade.  Especially,  the  PEMFC  offers  a  possible 
substitute  for  the  combustion  engine  owing  to  its  high  power  effi¬ 
ciency  and  low  level  of  pollution.  Although  most  of  the  studies  on 
fuel  cells  in  laboratories  are  focused  on  the  performance  of  a  sin¬ 
gle  cell  [  1—11  ],  they  are  very  helpful  for  us  to  understand  the  basic 
principle  of  the  fuel  cell  and  to  enhance  the  fuel  cell  performance. 
However,  for  a  large  fuel  cell  stack,  the  assembly  technique  is  one 
of  the  key  techniques  to  ensure  the  stack  having  a  long  lifetime  and 
high  reliability.  Previous  studies  have  shown  that  the  clamping  load 
(pressure)  in  the  assembly  of  a  single  cell  plays  an  important  role  in 
optimizing  the  performance  of  the  fuel  cell  [1-7].  However,  there 
has  been  no  detailed  study  for  the  assembly  of  a  large  fuel  cell  stack 
since  the  optimal  clamping  load  is  more  difficult  to  determine  than 
for  a  single  cell  due  to  the  complex  system  structure  and  working 
environment. 

The  clamping  load  of  a  large  fuel  cell  stack  affects  the  lifetime 
and  performance  of  the  stack  system  in  several  ways.  Too  large  a 
clamping  load  may  cause  some  components  in  the  stack  to  produce 
a  stress  high  enough  to  give  rise  to  plastic  deformation  even  cracks 
[12],  while  an  unreasonable  small  clamping  load  may  produce  a 
high  contact  electrical  resistance  at  the  interface  of  the  gas  diffu¬ 
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sion  layer  (GDL)  and  the  bipolar  plate  (BPP)  [1-3],  and  also  may 
cause  leakage  of  water  or  fuel  in  the  seal  interfaces  [13].  It  should 
also  be  pointed  out  that  even  though  a  reasonable  initial  clamping 
load  had  been  given  in  the  assembly  process  at  room  temperature, 
either  the  contact  electrical  resistance  or  the  seal  interface  pres¬ 
sure  or  the  structure  stress  may  be  out  of  the  allowed  region  of  the 
design  due  to  the  thermal  deformations  of  the  fuel  cell  components 
under  working  conditions.  To  overcome  those  problems  occurred 
during  the  stack  running,  the  stiffness  of  both  the  clamping  com¬ 
ponents  and  the  clamped  components  have  to  be  optimized.  For 
example,  in  the  bolt  assembling  (a  traditional  assembly  technique), 
both  the  stiffness  and  the  strength  of  the  bolts  should  be  first  opti¬ 
mized.  However,  the  following  problems  are  still  the  challenging 
techniques  in  the  structure  design  of  a  large  fuel  cell  stack:  (a)  how 
to  determine  the  clamping  load  of  the  stack;  (b)  how  to  design 
the  stiffness  and  strength  of  the  clamped  components;  (c)  how  to 
design  the  geometry  and  stiffness  of  the  sealants;  (d)  how  to  design 
the  bipolar  plates,  including  the  structural  shape  and  size,  etc. 

Strictly  speaking,  the  design  of  a  reasonable  clamping  load  and 
an  optimal  stiffness  of  the  related  components  for  a  large  fuel  cell 
stack  needs  a  detailed  finite  element  analysis  (FEA)  for  the  whole 
system.  Nowadays,  complex  FEA  can  be  obtained  by  virtue  of  highly 
developed  computer  technology  and  constantly  improved  software 
codes.  However,  it  is  difficult  to  form  an  ideal  FEA  model  of  the 
large  fuel  cell  stack  due  to  its  multi-scale  structural  characteris¬ 
tics.  For  example,  in  order  to  guarantee  the  reasonable  calculation 
accuracy,  for  the  membrane  electrode  assembly  (MEA),  the  size 
of  the  elements  should  be  on  the  order  of  micrometers  [1-3];  for 
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Nomenclature 

k  equivalent  stiffness  of  a  component  (N  mm-1 ) 

F  load  applied  on  the  component  (N) 

E  elastic  modulus  (N  mm-2) 

A  cross  section  area  (mm2) 

L  length  of  the  component  along  the  clamping  load 

direction  (mm) 

C  number  of  the  clamping  bolts 

N  number  of  the  flow  channels  on  one  side  of  the  bipo¬ 

lar  plate 

M  number  of  the  cells  of  a  large  stack 

ZF  clamping  deformation  (mm) 

ZT  free  thermal  deformation  (mm) 

AT  temperature  changes  (°C) 

Greek  letters 

a  coefficient  of  heat  expansion  (°C_1 ) 


the  interface  between  the  GDL  and  the  ribs  of  the  BPP,  the  con¬ 
tact  element  size  should  be  on  the  order  of  10  pan;  for  the  bipolar 
plate,  the  element  size  should  range  between  10  pan  and  100  pm; 
however,  for  the  end  plates  (EP)  and  the  assembling  bolts,  the  ele¬ 
ment  size  may  range  between  100  pm  and  1  cm.  Thus,  the  FEA  for 
a  large  fuel  cell  stack  is  a  multi-scale  nonlinear  contact  analysis 
coupled  with  thermal  stress  analysis,  a  really  hard  topic  in  the  com¬ 
putational  mechanics.  The  number  of  the  total  elements  may  be 
on  the  order  of  several  millions  even  larger.  No  report  has  been 
found  so  far  for  such  a  huge  computation  for  fuel  cell  structure 
analysis. 

In  this  paper,  we  present  a  high  efficient  assembly  technique 
for  the  large  PEMFC  stacks.  The  basic  idea  of  the  technique  is  first 
to  assume  the  stack  system  consisting  of  a  large  number  of  elas¬ 
tic  elements  (springs)  in  either  series  or  parallel  connections.  For 
each  spring,  the  stiffness  coefficient  can  be  obtained  using  material 
mechanics  analysis  for  the  simple  period  structures  or  using  FEA  for 
the  complex  structures.  Then  a  system  analysis  technique  is  used  to 
give  the  optimal  design  for  each  structure  component  in  the  stack. 

2.  Basic  model  descriptions 

Stiffness  of  a  solid  body  is  the  resistance  offered  by  an  elastic 
body  to  deformation.  Its  value  depends  on  the  material  property, 
the  structure  shape  and  size,  and  the  boundary  conditions  of  the 
solid  body  [14].  As  shown  in  Fig.  1,  for  an  elastic  bar  under  tension 
or  compression,  the  axial  stiffness  is: 


*  U 

F  F 


I  V  f  V  f '  fj  1 


Fig.  1.  Schematic  of  the  axial  equivalent  stiffness  of  an  elastic  element. 


where  k  is  usually  called  the  equivalent  stiffness  [15],  E  is  the  elastic 
modulus,  A  is  the  cross  section  area,  L  is  the  length  of  the  bar  along 
the  axial  direction.  As  for  a  complex  structure,  the  equivalent  stiff¬ 
ness  cannot  be  simply  obtained  by  Eq.  (1).  A  traditional  method  is 
to  simplify  the  structure  as  an  equivalent  model  consisting  of  elas¬ 
tic  elements  (usually  springs)  in  parallel  or  series  connections,  and 
then  the  deformation  of  the  structure  under  a  load  can  be  obtained 
through  this  equivalent  model.  In  this  paper,  the  most  we  care  about 
is  the  mechanical  behavior  of  the  fuel  cell  stack  along  the  direction 
of  the  clamping  load.  In  the  same  way,  the  whole  fuel  cell  stack  can 
be  divided  into  several  basic  components  whose  stiffness  can  be 
calculated  by  Eq.  (1)  according  to  the  load  bearing  characteristics 
of  the  structures  along  the  direction  of  the  clamping  load.  In  order 
to  make  the  model  building  process  more  comprehensive,  only  a 
single  cell  is  studied  in  the  first  step,  the  whole  stack  model  can  be 
obtained  via  cyclically  assembling  the  single  cells.  Thus  the  relation 
of  the  deformation  and  the  clamping  load  can  be  obtained  using  this 
equivalent  stiffness  model. 

2.1.  Basic  assumptions 

Proper  assumptions  are  needed  before  building  the  equivalent 
stiffness  model.  For  simplification  and  theory  analysis  strictness, 
the  following  assumptions  are  stated: 

(1)  The  end-plates  are  treated  as  rigid  body  and  have  an  enough 
high  mechanical  strength  [16].  Although  an  ideal  rigid  end-plate 
is  difficult  to  obtain  in  practice,  an  end-plate  with  a  suitable 
design  can  be  considered  as  rigid  one  approximately. 

(2)  Each  clamping  bolt  bears  the  same  clamping  load  and  gives  the 
same  elongation  during  assembly. 


Fig.  2.  Schematic  of  the  three  regions  of  the  bipolar  plate. 
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Table  1 

Stiffness  and  thermal  expansion  equations  of  the  basic  components. 
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Component 

Schematic 

Equivalent  Stiffness 

Thermal  Expansion 

BPP-base 

\F 

3 _ 

F  A 

^  BPP-base  ^  BPP-base 

A BPP-base 

a BPP-base  ^  BPP-base  ^ BPP-base 

BPP-rib 

fr-  g  1L  Atrr? 

L  ‘'Btr-nt  ,  _ 

'1 

E  BPP-rib  A  BPP-rib 

L BPP-rib 

&  BPP-rib^  T BPP-rib  ^  BPP-rib 

(3 )  The  bipolar  plate  has  N  parallel  flow  channels,  while  the  number 
of  the  bipolar  plate  ribs  is  (N  -  1 ). 

(4)  As  shown  in  Fig.  2,  the  bipolar  plate  is  divided  into  three  regions: 
the  BPP-sealant,  the  BPP-rib  and  the  BPP-base. 

(5)  MEA  is  considered  as  the  sandwich  structure:  GDL  +  PEM  (pro¬ 
ton  exchange  membrane)  +GDL  [17]. 

(6)  The  coolant  plate’s  stiffness  is  mixed  together  with  the  bipolar 
plate  as  a  single  bipolar  plate. 


(7)  It  is  assumed  that  there  are  no  gaps  at  the  connection  inter¬ 
faces  of  the  fuel  cell  stack  before  assembly,  i.e.,  all  the  MEA  and 
sealants  have  the  same  thickness. 

(8)  The  assembly  process  is  accomplished  at  room  temperature. 

In  addition,  in  the  following  analysis,  the  model  components 
listed  in  Table  1  are  called  the  basic  components,  whose  equiva¬ 
lent  stiffness  can  be  calculated  using  equations  given  in  Table  1. 


End  Plate 


Fig.  3.  Schematic  of  the  PEM  fuel  cell  stack  hardware. 


384 


P.  Lin  et  al.  /  Journal  of  Power  Sources  194  (2009 )  381  -390 


End  Plate 


K 

-71 BPP-base 


•  •••  h  •••• 

n  BPP -rib 


(N-l)  ribs 


•  •  •  • 


k  •  •  •  • 

nBPP-rib 


/i 


BPP-base 


h 


BPP-sealant 


u 

.  sealant 


Is 

>  <>  Z 


Tensile 
^  direction 

5  c  , 

>  —  y.K.,t 


BPP -sealant 


End  Plate 

~7 - 7 - 7 - 7 - 7 - 7 - 7 - 7 - 7 - 7 - 7 - 7 - 7^ 


Fig.  4.  Schematic  of  the  equivalent  stiffness  model  of  the  single  cell  (the  dimensions  of  the  springs  do  not  represent  the  magnitude  of  stiffness). 


The  components  consisting  of  the  basic  components  are  called  the 
assembled  components. 

2.2.  Equivalent  stiffness  model  of  a  single  PEMFC 

As  Fig.  3  depicts,  a  single  PEMFC  consists  of  three  types  of  com¬ 
ponents:  a  membrane  electrode  assembly,  two  bipolar  plates  and 
the  sealants  [17].  They  are  all  included  in  the  following  equivalent 
stiffness  model  of  the  single  cell. 

From  the  view  perpendicular  to  the  flow  channels,  we  have 
a  mechanical  equivalent  stiffness  model  of  the  single  cell  as 
shown  in  Fig.  4,  where  8  equals  the  thread  pitch  times  the 
pretightening  cycles  of  the  nuts,  and  the  arrow  Z  represents 
the  positive  direction  of  the  vectors  such  as  displacement  or 
load.  The  stiffness  of  each  spring  in  the  model  equals  the 
equivalent  stiffness  of  the  basic  component  or  the  assembled 
component. 

Based  on  assumption  (2),  we  can  make  further  simplification  for 
the  system  as  shown  in  Fig.  5.  This  model  is  a  pure  mathematical 
model  without  any  structural  characteristics,  the  guide  rails  will 
constrain  the  rotation  of  the  upper  end-plate  caused  by  the  unbal¬ 
anced  stiffness  distribution.  Any  way,  the  displacement  along  the 
clamping  direction  is  completely  equivalent  with  the  model  as  is 
shown  in  Fig.  4. 


Subsequently,  the  equivalent  stiffness  model  can  be  finally 
simplified  by  combining  the  stiffness  of  the  BPP-sealant  and  the 
sealant  (basic  components)  in  series  as  the  assembled  compo¬ 
nent:  external-cell,  and  combining  the  stiffness  of  the  BPP-base 
and  cell-core  in  series  as  the  assembled  component:  internal¬ 
cell.  After  those  simplifications,  the  system  of  one  single  cell  can 
be  simplified  as  what  is  shown  in  Fig.  6.  It  consists  of  two  end- 
plates,  C  clamping  bolts,  the  internal-cell  and  the  external-cell  (see 
Fig.  6). 

Based  on  the  model  shown  in  Fig.  6,  the  equivalent  stiffness  of  a 
PEM  single  cell  is  given: 


^single-cell  —  ^internal-cell  T  kexternal-cell  (2) 

where  k  is  the  equivalent  stiffness  corresponding  to  the  compo¬ 
nent  as  indicated  in  the  subscript  (similarly  hereinafter).  For  the 
internal-cell,  an  assembled  component,  the  equivalent  stiffness  can 
be  calculated  as: 


1 


'^internal-cell  — 


2/kBPP-base  +  l/^cell  -a 


/cBPP-base  fCcell-core 
^BPP-base  +  2/<ceii_cor 


(3) 


where  the  equivalent  stiffness  of  the  assembled  component  cell- 
core  (see  Fig.  5)  is: 
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Fig.  5.  Schematic  of  the  simplified  equivalent  stiffness  model  of  the  single  cell. 
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^cell-i 


(Ei=1  ^BPP-rib)  ^MEA 


C'C* + vkMEA  (E;':>(BPP-rib)  + 


2/<iv 


(4) 


where  the  superscript  i  represents  the  number  of  the  ribs.  The 
equivalent  stiffness  of  the  MEA  is: 

^GDL^PEM 


^MEA  = 


1 


(5) 


2  /  /<GDL  +  V  fcpEM  ^GDL  +  2  /<PEM 

In  the  same  way,  for  the  assembled  component  external-cell,  the 
equivalent  stiffness  can  be  calculated  as  follows: 

•  1  ^BPP-sealant  ^sealant 


,cexternal-cell 


2/^BPP-sealant  +  1  /^sealant  ^BPP-sealant  +  2/<seaiant 


(6) 


Equivalent  stiffness  of  the  basic  components  can  be  calculated  using 
the  equations  given  in  Table  1.  Based  on  the  conditions  of  com¬ 
patible  displacement  and  static  load  equilibrium  in  Z  direction 
(clamping  direction),  the  structural  load  equilibrium  equation  can 
be  written  as  follows: 

Cfcbolt^  +  ^single-cell)  +  ^single-cellZsingle-cell  =  0  (7) 

where  the  deformation  of  the  whole  single  cell  under  the  clamping 
load  is: 

_  C/<bolt  o 


ZF 

single-cell 


C^bolt  +  ^single-cell 


(8) 


where  Zfingle_cell  is  a  negative  number  indicating  the  structures  are 
under  compression  since  S  is  positive  for  bolt  stretching.  Then  the 
clamping  loads  applied  on  the  single  cell  are: 

^internal-cell  =  ^internal-cell^single-i 


i-cell 

^bolt  ^internal-cell 


(9) 


(10) 


C/<bolt  +  ^internal-cell  +  ^external-cell 

^external-cell  =  ^external-cell^single-cell 

_ _ Cfcbolt' /^external-cell _ ^ 

^/<bolt  +  ^internal-cell  +  ^external-cell 

The  total  clamping  load  applied  on  the  single  cell  is  the  sum  of 
internal-cell  and  ^external-cell .  which  equals  the  load  applied  on  the 
clamping  bolts  (see  Fig.  6): 

^clamping  —  ^internal-cell  T  ^external-cell  (H) 
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2.3.  Equivalent  damping  stiffness  model  of  a  PEMFC  stad< 

The  PEMFC  stack  is  comprised  of  numerous  cells  connected  via 
bipolar  plates  as  is  shown  in  Fig.  3.  For  this  cyclic  cell  characteristic, 
stiffness  of  the  whole  stack  can  be  obtained  by  combining  the 
single  PEMFC  equivalent  model  in  sequence.  Fig.  7  depicts  the 
equivalent  clamping  stiffness  model  of  a  PEMFC  stack.  Therefore, 
the  equivalent  stiffness  model  of  the  stack  can  be  simplified  as 
the  similar  model  shown  in  Fig.  6.  But  here  it  is  composed  of 
the  following  components:  two  end  plates,  C  clamping  bolts,  the 
internal-stack  and  the  external-stack  (see  Fig.  7).  Now  the  stack 
stiffness  can  be  written: 


^stack  —  ^internal-stack  +  ^external- stack 


(12) 


where  the  stiffness  of  the  internal-stack  and  the  external-stack  can 
be  calculated  as  follows: 


'^internal-stack  — 


^external-stack  — 


P-baSe)  +  E'!l(VCl-core) 

1 


E^1(l/fe®P-sea,ant)  +  E;:i(V^alant) 


M  ,  i  nM) 


(13) 


(14) 


In  Eqs.  (12)  and  (13)  the  superscript  i  represents  the  ith 
bipolar  plate  while  j  represents  the  jth  cell.  Stiffness  of  the 
cell-core  can  be  obtained  using  Eq.  (4).  Based  on  the  condi¬ 
tions  of  compatible  displacement  and  static  load  equilibrium, 
the  deformation  of  the  whole  stack  under  the  clamping  load 
is: 


ZF 

stack  ‘ 


Ckbolt 


Ckbolt  +  kstack 


(15) 
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Also  the  clamping  loads  applied  on  the  stack  are: 

^internal-stack  =  ^internal-stack^stack 

_  ^bolt  ^internal-stack  o 


^bolt  T 


internal-stack 


+  k 


external-stack 


(16) 


^external-stack  —  ^external-stack^stack 

^bolt  ^external-stack 


Ckbolt  + 


internal-stack  +  kexternal- stack 


(17) 


The  total  clamping  load  applied  on  the  stack  is  the  sum  of 
^internal-stack  ^external-stack*  which  equals  the  load  applied  on  the 
clamping  bolts: 


^clamping  —  ^internal-stack  +  ^external-stack  (18) 

2.4.  Amendments  to  the  equivalent  stiffness  model  of  PEMFC  stack 

What  is  discussed  above  is  based  on  the  assumption  that  there  is 
no  gap  before  assembling  at  either  the  interface  between  the  bipolar 
plate  and  the  MEA  or  the  interface  between  the  bipolar  plate  and 
the  sealant.  This  is  only  an  ideal  case.  In  practice,  there  are  usually 
small  gaps  before  assembling  at  one  of  the  interfaces  mentioned 
above.  We  will  discuss  this  general  case  in  Section  2.4.1. 

The  clamping  process  is  usually  carried  out  at  room  temper¬ 
ature.  However,  the  modern  PEM  fuel  cell  stack  usually  operates 
at  temperature  below  80  °C  [18].  All  the  components  will  produce 
thermal  stress  and  deformation  under  such  a  temperature.  We  will 
also  discuss  the  temperature  effect  in  Section  2.4.2. 

2.4.1  Unequal  thickness  of  the  MEA  and  the  sealant  before 
assembly 

In  the  case  of  unequal  thickness  of  the  MEA  and  the  sealant 
before  assembly,  the  bipolar  plate  will  first  contact  with  either  the 
sealant  or  the  GDL.  It  is  assumed  that  all  the  cells  have  the  same 
gap,  Zfnitial,  defined  as  below: 

^initial  =  ^MEA  —  ^sealant  (19) 

where  LMEA  and  Lsealant are  the  thickness  of  the  MEA  and  the  sealant, 
respectively.  When  ZFnitial  is  positive,  the  gap  occurs  at  the  interface 
of  the  sealant  and  the  bipolar  plate;  but  when  Zfiitial  is  negative,  the 
gap  occurs  at  the  interface  of  the  MEA  and  the  bipolar  plate. 

Here  the  total  pretightening  displacement  of  the  clamping  bolts 
is: 


^total  —  ^initial  +  ^fasten 


(20) 


where  5initial  is  the  pretightening  displacement  occurring  before  the 
gap  vanishes  at  the  initial  stage  of  the  assembly,  5fasten  is  the  design 
displacement  variable  of  the  clamping  bolts  after  gaps  vanish. 

In  the  case  of  Emea  >LSealant,  we  have  (see  Fig.  8a): 


8 


initial  — 


^^initial^^bolt  ^internal-stack) 

Ckbolt 


(21) 


where  M  is  number  of  the  cells  of  a  large  stack  and  C  is  the  number 
of  the  bolts.  With  Eq.  (14),  the  total  deformation  of  whole  stack, 
ZtFotai*  can  be  obtained: 


Ztota>  =  —(^initial)  +  4ack  =  -(^initial) 


Ck, 


bolt 


Ckbolt  +  ^stack 


^fasten 


(22) 


where  Zjtack  denotes  the  whole  stack  deformation  after  the  gap 
vanishes.  Here  the  clamping  load  applied  on  the  stack  structure  is: 


(a) 


(b) 


BPP 

7F 

" initial 

/■ 

MEA 

sealant 

BPP 

BPP 

7f 

" initial 

MEA 

sealant 

BPP 

Fig.  8.  Schematic  of  unequal  thickness  of  the  MEA  and  sealant  before  assembly:  (a) 
a  gap  exists  at  the  interface  of  the  sealant  and  BBP;  (b)  a  gap  exists  at  the  interface 
of  the  MEA  and  BBP. 


^internal-stack  —  ^internal-stackZ^otai 


—  ^internal- 


stack 


(MZfniti 


+  ' 


Ckbolt 


initial  '  C/<bolt  +  /<stack 


r  _  VF  _  ^bolt ^external-stack 

r  external-stack  —  '^external-stack^stack  —  rb  i  b 

Ckbolt+kstack 

In  the  case  of  LMEA  > lseaiant* we  have  (see  Fig.  8b): 

c  ^^initial^^bolt  +  ^external-stack) 

^initial  = 


Ckbolt 

and  the  total  deformation  of  the  whole  stack  is: 
Ckboit 


Z'  ,  =  MZl 


^fasten 


■"total  initial  rb,  ,  \  b 

CK-bolt  T-  '<-stack 

Also  the  clamping  loads  applied  on  the  stack  are: 


^internal-stack  —  ^internal-stackZ^^  — 


^bolt  ^internal-stack 
^bolt  +  ^stack 


°fasten  J 

(23) 
^fasten  (24) 

(25) 

(26) 

^fasten  (27) 
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^external-stack  —  ^external-stack-^total 

=  ^external-stack  Ck^  +  k^' 


(-^bolt 


'fasten 


,)  (28) 


2.4.2.  Effect  of  temperature 

Temperature  variation  can  cause  thermal  deformation.  For  the 
structure  shown  in  Fig.  1,  when  the  temperature  is  changed  by  AT, 
the  free  axial  thermal  deformation  will  be: 


A  L  =  aATL 


(29) 


If  this  deformation  is  partially  or  fully  constrained,  the  thermal 
stress  will  be  generated  inside  the  fuel  cell  structures.  As  shown 
in  Fig.  9,  the  thermal  deformation  of  the  whole  fuel  cell  stack  can 
be  obtained  based  on  the  compatible  displacement  and  static  load 
equilibrium  conditions: 


C/<bolt^b0it  +  ^internal-stack^temal-stack 


ZT 

stack  ' 


+  /<( 


extemal-stack^Jxtemal-stack 


C^bolt  +  ^internal-stack  +  ^external-stack 


(30) 


According  to  Eqs.  (15)  and  (16)  the  combined  loads  applied  on  the 
stack  are: 


^internal-stack  —  ^internal-stack-^. 


stack 


+  k\ 


internal- 


stack  (  Zj 


internal-stack 


+  ZJtack) 


(31) 
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Fig.  9.  Schematic  of  the  equivalent  stiffness  model  of  the  PEM  stack  when  consid¬ 
ering  the  thermal  effect. 


^external-stack  —  ^external-stack-^stack 

+  ^external-stack  ( — ^external-stack  +  Zstacl<) 

where  the  free  thermal  deformation  of  the  internal-stack  and  the 
external-stack  can  be  given  using  the  equivalent  clamping  stiffness 


(32) 


model  shown  in  Fig.  7: 
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Fig.  10.  Dimensions  of  the  single  cell  model. 
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(34) 


In  addition,  the  thermal  deformations  of  the  basic  components  can 
be  calculated  referring  to  Table  1.  Therefore  the  thermal  defor¬ 
mations  of  the  assembled  components  are  obtained  based  on  the 
thermal  deformations  of  the  basic  components  (see  Fig.  4). 


2.5.  Determination  of  the  assembly  parameters  using  the 
equivalent  stiffness  model 

Systematical  use  of  the  equivalent  stiffness  model  can  form  a 
high  efficient  assembly  technique  for  fuel  cell  stacks. 

(1)  After  obtaining  the  basic  information  such  as  component 
dimensions  and  material  properties,  the  equivalent  clamping 
stiffness  of  the  basic  components  and  assembled  components 
can  be  calculated  using  the  method  mentioned  above. 

(2)  Based  on  the  obtained  stiffness,  the  clamping  load  of  the  stack 
and  the  contact  pressures  on  the  important  interfaces  are  cal¬ 
culated.  Thus  the  internal  stress  of  the  components  can  be 
calculated  and  the  contact  electrical  resistance  can  be  obtained 
using  the  equations  given  in  refs.  [1-3]. 

(3 )  The  range  of  the  clamping  load  can  be  determined.  For  the  lower 
bound,  the  clamping  load  should  be  large  enough  to  provide  rea¬ 
sonable  contact  pressures  for  the  contact  electrical  resistance 
and  the  sealing  pressure.  For  the  upper  bound,  the  clamping 
load  should  not  be  so  large  that  the  internal  stress  of  the  basic 
components  reaches  the  yield  stress.  Finally  the  optimal  clamp¬ 
ing  load  can  be  determined  and  other  clamping  parameters 
such  as  the  stiffness  and  strength  of  the  clamping  bolts  can  be 
designed  in  the  same  way. 


3.  Comparisons  of  the  model  analysis  and  FEA  analysis 

In  order  to  verify  the  equivalent  stiffness  model,  comparisons 
of  the  model  analysis  and  the  accurate  analysis  using  FEA  have 
been  carried  out  based  on  the  PEM  single  cell  described  in  refs. 
[5,19-21  ].  The  commercial  FEA  code  of  MSC.  Marc  was  used  in  this 
study. 

3.1.  FEA  model  definition 

3.1.1.  Geometry 

Based  on  the  two-dimensional  (2  D)  FEA  model  proposed  in  refs. 
[5,19-21  ],  a  three-dimensional  (3  D)  FEA  model  is  established.  It  is 
1  /8  of  a  single  cell  due  to  the  structural  symmetry.  The  present  3 
D  model  consists  of  the  bipolar  plate,  the  GDL,  the  PEM  and  the 
sealant,  while  the  end-plates  and  the  clamping  bolts  are  not  con¬ 
sidered  in  the  numerical  analysis.  Dimensions  of  the  PEM  single  cell 
can  be  obtained  from  Fig.  10,  other  parameters  not  defined  in  the 
figure  are:  the  number  of  the  ribs  (35),  the  thickness  of  the  PEM 
(0.05  mm),  the  thickness  of  the  GDL  (0.275  mm),  and  the  thickness 
of  the  sealant  (0.3  mm).  Here  the  MEA  and  sealant  have  the  same 
thickness. 

3.1.2.  Mesh  and  mechanical  properties 

Mapped  meshing  is  adopted  in  order  to  ensure  proper  ele¬ 
ment  connectivity  and  reasonable  aspect  ratio.  The  complete 
mesh  includes  240,211  nodes  and  195,516  elements  (see  Fig.  11a). 
Full  integration  hexahedron  elements  (element  7)  have  been 
used.  Mechanical  properties  of  each  component  are  listed  in 
Table  2.  All  the  components  are  assumed  to  have  linear  elastic 
behavior. 
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Fig.  11.  Three-dimensional  contact  model:  (a)  mesh  and  (b)  contour  plot  of  the 
deformation  along  the  clamping  direction  (z-axis). 
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Fig.  12.  Relations  between  the  load  and  the  displacements  of  the  single  cell  structures. 


3.1.3.  Boundary  conditions 

Symmetry  constraints  are  set  on  the  symmetrical  planes  in  order 
to  prevent  free  movements.  The  efficient  contact  algorithm  of  MSC. 
Marc  is  used  to  compute  the  interactions  between  the  components. 
In  addition,  the  end-plate  is  replaced  by  a  rigid  surface  controlled 
by  the  control  node  [22].  The  clamping  load  increasing  linearly  from 
0  to  1865  N  per  bolt  [5]  is  directly  applied  on  the  control  node  to 
give  a  uniform  pressure  on  the  bipolar  plate. 

3.2.  Simulation  results  and  comparisons 

Fig.  lib  shows  that  under  the  uniform  pressure,  all  the  com¬ 
ponents  deform  nearly  uniformly  along  the  clamping  direction 
according  to  the  well-layered  contour  plot.  The  contour  lines  only 
fluctuate  slightly  at  structural  junctions  especially  for  the  bipolar 
plate  consisting  of  BPP-base,  BPP-rib  and  BPP-sealant.  This  phe¬ 
nomenon  strongly  supports  the  initial  assumption  (4).  Using  the 
control  node,  displacement  of  the  single  cell  structure  can  be 
obtained.  In  the  same  way,  a  load  deformation  diagram  of  the 
1/8  single  cell  structure  is  given  in  Fig.  12.  It  shows  that  both  the 
single  cell  and  its  components  deform  almost  linearly  with  the  dis¬ 
placement.  Therefore  the  corresponding  equivalent  stiffness  can  be 
obtained  by  calculating  the  slopes  of  the  curves. 

Meanwhile,  the  corresponding  equivalent  stiffness  of  the  com¬ 
ponents  can  be  calculated  using  the  equations  mentioned  above 
and  those  given  in  Table  1.  In  addition,  the  two  layers  of  sealants  of 
the  single  cell  can  be  combined  in  series  as  a  whole  in  the  model 
analysis  (MA).  Thus,  for  the  same  structure,  the  MA  is  compared 
with  the  FEA  in  Table  3.  The  relative  error  using  the  model  analysis 
is 

abs(MA-FEA) 

Error  =  — - — — — - -  x  100%  (35) 

FEA 


Table  2 

Material  properties  of  the  single  cell  hardware. 


Table  3 

Result  comparisons  of  the  MA  and  the  FEA. 


Components 

Equivalent  stiffness 
(xlO6  N  mm'1)  MA 

Equivalent  stiffness 
( x  106  N  mm-1 )  FEA 

Relative 

error 

BPP-base  (1/4)* 

0.523324 

0.523325 

<1 

BPP-rib  (1/2) 

0.17395 

0.17395 

<1 

BPP-sealant  (1/4) 

0.466406 

0.46641 

<1 

Sealant  (1/4) 

10.945 

10.945 

<1 

GDL  (1/4) 

22.45536 

22.45534 

<1 

PEM  (1/8) 

4.69316 

4.69318 

<1 

Single  cell  (1/8) 

0.84783 

0.83961 

<3 

Values  in  the  brackets  in  column  1  indicate  the  symmetry  of  the  basic  compo¬ 
nents  and  the  single  cell  structure. 


Table  3  shows  that  the  results  calculated  by  the  two  methods  are 
very  close.  The  relative  error  for  the  basic  component  is  less  than  1  %, 
and  even  for  the  whole  single  cell  structure  the  relative  error  is  no 
more  than  3%.  In  addition,  it  costs  nearly  10  h  for  the  FEA  to  obtain 
the  results  on  a  high  performance  hardware  platform  with  4  cores 
CPU  and  4  GB  memory;  while  it  costs  about  15  min  by  hand  cal¬ 
culation  and  it  may  need  only  a  few  seconds  using  the  computer 
programming  based  on  this  model  analysis.  All  the  calculations 
mentioned  above  show  that  the  equivalent  stiffness  model  has  a 
highly  efficiency  and  gives  a  very  good  prediction  accuracy. 

4.  Discussions  and  conclusions 

The  clamping  load  applied  on  the  fuel  cell  stack  should  be  fit  for 
the  corresponding  structures  since  the  mechanical  behavior  of  the 
fuel  cell  stack  has  a  profound  affect  on  the  fuel  cell  performance.  Too 
high  a  clamping  load  will  not  only  reduce  the  permeability  of  the 
GDL  and  then  reduce  the  power  density  of  the  fuel  cell  stack,  but 
also  may  cause  some  components  in  the  stack  to  produce  a  high 
stress  as  a  result  that  a  plastic  deformation  even  crack  may  take 
place;  Flowever,  an  unreasonably  low  clamping  load  may  cause  a 
high  contact  electrical  resistance  at  the  interface  of  the  GDL  and 
the  BPP,  and  a  leakage  of  either  water  or  fuel  in  the  seal  interfaces 
[1-3].  Moreover,  even  though  the  average  stress  in  a  component  is 
lower  than  the  yield  stress,  a  fatigue  failure  may  occur  during  the 
running  time  of  the  stack  if  the  stress  variation  range  due  to  either 
the  temperature  or  other  working  conditions  is  considerably  large. 

In  this  paper,  a  large  fuel  cell  stack  has  been  simplified  to  an 
equivalent  stiffness  model  consisting  of  numerous  elastic  elements 
(springs)  in  parallel  or  series  connections.  Starting  with  a  single 
PEM  fuel  cell,  an  equivalent  stiffness  model  of  the  whole  fuel  cell 
stack  is  developed.  Then  the  effects  of  the  temperature  and  struc¬ 
ture  parameters  on  the  internal  stress,  contact  resistance  and  the 
sealing  pressure  of  the  components  are  discussed.  Finally,  a  high 
efficient  assembly  technique  for  a  large  PEMFC  stack  using  the 
equivalent  stiffness  model  has  been  proposed.  With  this  technique, 
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a  reasonable  clamping  load  and  optimal  stiffness  of  the  related 
components  for  a  large  fuel  cell  stack  can  be  obtained. 

Comparisons  of  the  model  analysis  and  the  accurate  analysis 
using  FEA  have  been  carried  out  based  on  a  typical  PEM  single  cell. 
A  numerical  analysis  using  a  three-dimensional  nonlinear  contact 
model  shows  that  the  prediction  error  of  the  stiffness  for  each  com¬ 
ponent  of  a  single  fuel  cell  using  the  present  model  is  less  than  1%, 
and  even  for  the  whole  single  cell  structure  the  relative  error  is  no 
more  than  3%. 

In  the  present  part  of  the  study  for  the  assembly  technique  of 
a  large  fuel  cell  stack,  we  used  a  rigid  end-plate  assumption.  This 
may  causes  some  error  for  the  structure  design  and  the  prediction  of 
the  internal  stress  and  the  interface  contact  pressure  of  the  fuel  cell 
components,  especially  for  the  cells  near  the  end-plates.  However, 
if  the  end-plate  is  so  specially  designed  that  the  interfaces  of  the 
components  inside  the  stack  can  be  kept  as  an  approximate  plane, 
the  rigid  end-plate  assumption  is  still  valid.  The  work  how  to  modify 
the  present  assembly  model  for  an  elastic  end-plate  and  how  to  use 
the  present  model  to  design  the  structures  for  a  large  fuel  cell  stack 
is  being  done  in  the  present  research  group  and  will  be  published 
in  additional  paper. 
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